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Abstract 

An introduction to the defects which dominate the physics of superfluid 
He 4 films, of superconducting slabs and of crystalline and hexatic membranes 
is given. We first review point vortices in two-dimensional neutral super- 
fluids and discuss the unusual screening which arises when the bosons are 
charged, as in superconducting films. Dislocation and disclination defects in 
crystalline membranes are discussed from a similar point of view. There is 
little or no screening in "monolayer" crystals, which are strongly constrained 
to lie in a flat two-dimensional plane. A strong nonlinear screening effect 
arises, however, in 2d membranes allowed to buckle into the third dimension. 
This screening drastically lowers dislocation and disclination energies, and 
forces crystalline membranes to melt at any finite temperature. We point out 
that buckled 5- and 7-fold disclinations in hexatic membranes have in gen- 
eral different logarithmically divergent energies. A similar asymmetry exists 
in the energies of 5- and 7- fold defects in liquid membranes. This difference 
determines the sign of the Gaussian bending rigidity, and has important conse- 
quences in membranes which can change their topology or with free boundary 
conditions. 



* Lectures given at the NATO ASI on Fluctuating Geometries in Statistical 
Mechanics and Field Theory, Les Houches, Summer 1994. 
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I. INTRODUCTION 



Phase transitions and spontaneous symmetry breakings abound in modern condensed 
matter physics. Of particular interest are the low energy excitations associated with con- 
tinuous broken symmetries, like the broken translational and orientational invariance of a 
crystal, the broken rotational symmetry of ferromagnet or the broken gauge symmetry of 
superfluids and superconductors. By imposing a slow spatial variation of the continuous 
symmetry operation on the ground state, one generates an important class of modes which 
can be easily excited near T = 0. In ferro magnets for example, these Goldstone modes are 
called spin waves, or, in quantized form, "magnons." In crystalline solids, the translational 
broken symmetryl leads to new shear phonon modes. Unlike a nematic liquid crystal, the ad- 
ditional broken rotational symmetry of a solid does not lead to any extra low energy modes: 
Rotational modes are coupled to phonons in a way reminiscent of the Higgs mechanism and 
are thus rendered massive — see Ref. || below and references therein. In superfluid helium, 
the broken symmetry selects a particular phase for the condensate order parameter, and 
leads to second sound. The physics of superconductors is complicated by the charge of the 
Cooper pairs and by the coupling of electron motion to the underlying ionic lattice. Low 
energy phase-degrees of freedom are still important, although now it is the vector potential 
which acquires a mass via the usual Higgs mechanism. 

The excitations discussed above are all continuously connected to the ground state. How- 
ever, point, line and wall defects in the order parameter texture, which are not continuously 
connected to the ground state, may also be important. Dislocation lines mediate plastic 
flow, and point-like vacancies and interstitials allow rapid particle diffusion in crystalline 
solids. Vortex rings control the decay of supercurrents in superfluids and superconductors. 
Most laboratory ferromagnets are subject to weak crystal fields, which break rotational in- 
variance. In this case, the width of hysterisis loops below the Curie temperature is controlled 
by the motion of bloch walls. Under some circumstances, such topological defects also play 
an important role in the phase transition out of the broken symmetry state. 

In this review, we discuss the defects which arise in superfluids, superconductors and 
membranes. Although the emphasis is how to incorporate defects into simple Landau- 
Ginzburg field theories, we also touch upon their role in the statistical mechanics. We 
concentrate on two-dimensional systems, so that the defects are typically point-like. 

The classic example of vortices and the Kosterlitz-Thouless transition in superfluid he- 
lium films is discussed first, together with universal jump in the superfluid density. Experi- 
mental confirmation of this prediction was an important test of the idea of defect-mediated 
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phase transitions. Screening of defect energies in superconducting films is then treated, with 
particular attention to the spreading of magnetic field lines in the three-dimensional volume 
outside the sample. When this spreading is ignored, screening efficiently cuts off the loga- 
rithmic divergence in the vortex energy. The finite energy vortices in this "naive" theory 
then proliferate for entropic reasons at any finite temperature, thus destroying supercon- 
ductivity. In realistic samples, however, magnetic field line spreading (which arises because 
supercurrents are confined to a two-dimensional plane) renders screening less efficient, and 
a Kosterlitz-Thouless transition from a finite temperature superconducting state becomes 
possible. 

The remainder of the paper contrasts the behavior of defects in membranes and monolay- 
ers. By "monolayers," we mean crystalline films strongly confined to a plane. Dislocations 
and disclinations control the statistical mechanics of monolayers in a way similar to vortices 
in neutral superfluids. Successive dislocation and disclination unbinding transitions lead 
from crystalline to hexatic to liquid monolayer phases. In contrast to monolayers, "mem- 
branes" have a relatively low energy cost to buckle out of the two-dimensional plane, thus 
escaping into the third dimension. As a result, defects can be screened by out-of-plane dis- 
placements much like vortices in the "naive" model of superconducting films. Dislocations, 
for example, now have a finite energy. Their entropic proliferation causes crystals to be 
replaced by hexatics as the inevitable low temperature phase of membranes. We also discuss 
the buckling of 5- and 7-fold disclinations, of vacancies, interstitial and impurity atoms, and 
of grain boundaries. 

We conclude by emphasizing the fundamental asymmetry between 5- and 7-fold disclina- 
tions in membranes. There is an exact symmetry between plus and minus vortices in super- 
fluids and between dislocations with equal and opposite burgers vectors in two-dimensional 
crystals, which forces the energies of the defect and anti-defects to be equal. No such sym- 
metry, however, insures that the core energies of plus and minus disclinations are identical. 
These disclinations represent points of local 5- and 7-fold coordination, where the number of 
neighboring particles is determined by the Dirichlet construction. This situation is similar 
to vacancies and interstitials in crystalline solids: Vacancies typically have a lower energy 
than interstitials, even though these excitations are "anti-defects" of each other. 

Disclination asymmetry is of little consequence in flat monolayers, where the concentra- 
tions of 5- and 7-fold defects must be equal for topological reasons. (More precisely, Euler's 
theorem ensures that the average coordination is six in a sufficiently large system.) In hex- 
atic monolayers, moreover, the logarithmic divergences in the energies of isolated plus and 
minus disclinations have equal coefficients. The situation is different, however, in membranes 
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which are free to change their topology or with free boundary conditions which allow them to 
curl up near the edges. When the bending rigidity which controls out-of-plane fluctuations 
is small, buckling in hexatic membranes leads to different coefficients in the logarithmically 
diverging disclination energies. In liquid membranes, moreover, different core energies of 
plus and minus disclinations will inevitably bias the system towards a net positive or nega- 
tive average Gaussian curvature. The energy asymmetry between 5- and 7-fold disclinations 
in the liquid determines the sign of the Gaussian rigidity, which acts as chemical potential 
controlling their population difference. Returning to the comparison with a superconductor, 
the Gaussian rigidity acts like an external magnetic field, which would bias superconducting 
films toward a net excess of plus or minus vortices. In membranes, the sign of the Gaussian 
bending rigidity determines whether vesicle phases (with positive net Gaussian curvature) 
or "plumbers nightmare" phases (with overall negative Gaussian curvature) are preferred. 

Theories and simulations of hexatic and liquid membranes usually avoid this issue by 
imposing a definite spherical or torroidal topology, which automatically fixes the difference 
in the number of 5- and 7-fold defects. Such constraints are clearly artificial, however, 
for lipid bilayer membranes which can change their topology at will, and may even have 
free edges at some intermediate stage of a topological change. The analogous constraint 
for vacancies and interstitials in solids would be imposed by periodic boundary conditions: 
these defects can then only be created in pairs. In real solids, with free boundary conditions, 
vacancies and interstitials can enter at the surface to produce the concentrations dictated 
by the microscopic energetics. The simplest probe of disclination bias would be to simulate 
liquid membranes with free edges, and see which defect begins to dominate as the bending 
rigidity is reduced from large values. See section III.D for details. 

II. TWO-DIMENSIONAL SUPERFLUIDS AND SUPERCONDUCTORS 

A detailed and successful theory of two-dimensional superfluidity was constructed in 
the 1970s. The theory explains the observed properties of thin films of superfluid He 4 on 
various substrates and the physics of thin superconducting films. The basic physical idea 
of defect unbinding was proposed by Kosterlitz and Thouless JI|], by Berezinski [0, and, in 
a different physical context, by deGennes ||. The first detailed calculations were carried 
out by Kosterlitz M, who exploited a pioneering renormalization group method developed 
by Anderson and Yuval for the Kondo problem [[J. Because extensive reviews are available 
elsewhere we limit ourselves to a brief sketch of the models and the basic results. 

We shall see how vortices control the basic physics and discuss their statistical mechanics. 
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We discuss only equilibrium properties, and ignore the important area of superfluid and 
superconducting dynamics [P|,pHf|. 



Oscillating Substrate 



He 4 



(a) 



Third Sound 



(b) 

FIG. 1. (a) Schematic of an oscillating substrate experiment, which 
measures the superfluid density, (b) Propagation of a third sound excitation. 
Measuring the dispersion relation also allows one to extract p s (T). 



A. Superfluid Helium Films: Experimental Facts 

Superfluidity arises with decreasing temperature in helium films in a fundamentally dif- 
ferent way than in bulk materials. Two basic experiments which reveal this difference are 
illustrated in Fig. 1. On an oscillating substrate, a finite fraction of the He 4 film decouples 
from the underlying motion below a critical temperature T c . The film appears lighter than 
it should be, and measuring this deficit allows one to extract the temperature-dependent 
superfluid density p s {T) ||11|| . In a related experiment, third-sound excitations similar to 
capillary waves propagate below T c with negligible damping through films as thin as a frac- 
tion of an Angstrom |T2|] . Such propagating waves would be impossible in a classical fluid, 



due to viscous damping. The primary restoring force is the Van der Waals attraction to the 
substrate (rather than gravity or surface tension), and again one can extract p s {T) provided 
the temperature is not too close to T c , where the waves become strongly damped. 
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FIG. 2. Jump discontinuity in the superfluid density p s {T) as a function 
of temperature for a fixed thickness helium film. Also shown is the superfluid 
specific heat at constant pressure as a function of temperature. There is 
only an essential singularity, C P (T) m c\ + C2exp[— c$/\T — T c | 1//2 ], near the 
transition, where the Cj are constants. 

As illustrated in Fig. 2, the superfluid density in films actually jumps discontinuously 
to zero at T c , in contrast to the continuous vanishing observed in bulk superfluid He 4 . The 
phase transition is nevertheless not first order, as becomes evident from an examination 
of the specific heat (see Fig. 3). There is only a rounded maximum above the transition 
temperature, in contrast to the delta function peak expected for a true first-order transition. 
The detailed theory shows that the specific heat contains only an unobservable essential 
singularity at the transition itself |14],T3| . 



B. Theoretical Background 

Understanding the equilibrium behavior of superfluid helium films at temperature T 
requires calculation of quantities such as the quantum partition function, 

Z = Tr'{e-^ T ] , (2.1) 

where H. is the Hamiltonian operator, 



-h 2 v 

2m 



£ + (2.2) 



j=l i>j 



and Tr' means a sum over all symmetrized energy eigenfunctions, reflecting the bosonic 
nature of a He 4 atom. Here V(r) is the boson pair potential, and we have set Boltzmann's 



6 



constant ks = 1- Explicit computation of Z for a fixed large number of particles N by 
finding all symmetric energy eigenfunctions (or even just the ground-state wave function 
|15| ) is difficult and it is useful to bypass this problem by changing to a second quantized 



representation of the many body physics. In second quantized notation, the Hamiltonian 
reads 

H = Jd 3 r $+{r) V 2 ^) ^(f) 

+ \i ^ r I ^ - r'\)$(r')ip(r) , (2.3) 

where ?/> + (r) and ijj(r) are boson creation and destruction operators obeying the commutation 
relations, 

[^(f),^ + (f')] = 5(f-r '), (2.4) 

with all other commutators zero. Within the second quantization formalism, it is convenient 
to calculate thermodynamic quantities via the grand canonical partition function, 



Z 



Tr > L-{h-»n)/t\ (2 5) 



where /i is the chemical potential and the number operator is 

N= f d 3 r ^ + (f)^(r) . (2.6) 



As pointed out by Penrose and Onsager [0, Bose condensation in bulk He 4 is associated 



with off-diagonal long-range order in the correlation function, 

G(r) = ^ + (r>(C))>, (2.7) 

where (Q) means Tr'{Q exp[— (H — /jN)/T]}/Z gr . Specifically, one has 

;imG(f) = (^+(f)) $(3)) 

^ . (2.8) 

A simplified description of superfluidity in terms of a Landau theory results from a coarse 
graining procedure. In the case of helium films, we start by defining a c-number field 

i>(f) = (4>) m , (2.9) 

where the brackets include a spatial average over a volume Q(f) centered on r which is large 
compared to the spacing between helium atoms but small compared to the overall system 
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size. The coarse-grained fields ip*(r) and ip(r) (which are nonzero when the He 4 liquid is 
superfluid in view of Eq. ( [2.8|) ) behave like complex numbers (instead of operators) for large 
enough averaging volumes. In helium films, we also take the averaging size to be large 
compared to the film thickness. The long wavelength spatial configurations of the system 
are now specified by complexions of ip(r), where r is a £wo-dimensional spatial variable, and 
the partition function is a functional integral, 

Z gr = J 2^(f)exp[-F/T] . (2.10) 

The free energy F includes entropic contributions due to the coarse-graining procedure. The 
Landau expansion of F/T in the order parameter ip(r) has the same form as Tt — pN, with 
ip(r) and ?/> + (r) replaced by the c-number fields ip(r) and tp*(r), 



F 
T 



d 2 r 



l -A\V^\ 2 + l -a\^\ 2 + b\^ + 



(2.11) 



where, as usual, a changes sign at the mean-field transition temperature T c °, a = a'(T — T c °). 

Because the long wavelength fluctuations in a helium film are two-dimensional, the actual 
temperature T c where true long-range order develops is suppressed well below T c °. For 
T pa T c , a <C and the polynomial part of F has a deep minimum, even though thermal 
fluctuations may nevertheless suppress genuine diagonal long-range order at the largest 
length scales. To a first approximation, we can then neglect fluctuations in the amplitude of 
ip(r), and set 

</>(f) pa ifj e ie ^, (2.12) 



where ipo pa ^J—a/Ab and 6{r) is a slowly varying phase variable. The free energy becomes 

I = const . + l Ko j |w(r)| 2 d\ (2.13) 

where K = AipQ. The physical interpretation of K Q becomes clear if we recall that the 
superfluid velocity v s is given by [TT 



v s (r) = — W(f), (2.14) 
m 

and write the contribution to F of the suprfluid kinetic energy in terms of the superfluid 
density p° s as 

F = const. + ^p° s J d 2 r\v s (r)\ 2 . (2.15) 
Evidently, K is related to the superfluid density in this "phase only" approximation by 
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K ° = St p ° • (2 - 16) 

The peculiar nature of superfluid order in two dimensions becomes evident if we now 
evaluate G(r), by integrating freely over the phase field when evaluating averages weighted 
by e~ F / T , 



G( r ) = (mno)) 



A exp 



(2.17) 



where the last line is a general property of Gaussian fluctuations. Upon introducing Fourier 
variables, the remaining average is easily evaluated using the equipartition theorem, 

{m - mn = A/|^ (1 _, f -, 

r 5»7^ ln(r/o ° ) - (2 - 18) 



where a is a microscopic cutoff. We thus find that G(r) decays algebraically to zero fig] , 

G(r) ~ , (2.19) 

with 

r?(T) = l/27tK (T) . (2.20) 

Thus there is only quasi-long range off-diagonal long-range order at any finite temperature 
in two dimensions. Nevertheless, this algebraic order is enough to distinguish the low- 
temperature superfluid from a distinct high temperature normal liquid film in which G(r) 
decays exponentially, 

G(r) ~ ex P [-r/£(T)] , (2.21) 

which defines a superfluid coherence length £(T). Note that the low-temperature phase 
still has a nonzero superfluidity density, similar to bulk superfluids. In contrast to bulk 
superfluids, however, the low-temperature phase is characterized by a continuously varying 
critical exponent rj(T). 
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(a) 





V(*,0)|/y 
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(b) 

FIG. 3. (a) Lines of constant phase associated with a vortex in a two- 
dimensional superfluid. Dashed circles are the streamlines for the associated 
supercurrents. (b) Variation of the order parameter magnitude near a vortex 
core located at the origin. 

The "phase only" description used above neglects amplitude fluctuations, and ignores 
the periodic nature of the phase variable. We can account for both effects by introducing a 
discrete set of vortex singularities in if>(r). As illustrated in Fig. 3a, vortices represent zeroes 
of ip(r), where the phase is undetermined. Near a vortex core, the assumption that the order 
parameter amplitude ^(r)! is a constant breaks down. Fig. 3b shows how ip(f) varies in the 
vicinity of a vortex ||19|| . ip(r) rises from zero to its bulk value over a characteristic length 



£0; related to the coefficients of the first two terms of Eq. ( |2.11| ), 

£ = v^Th ■ 



(2.22) 



To describe a radially symmetric net phase change of ±27r near a vortex centered at the 
origin, we set 
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x,y) = ±taiT 1 (y/x) . (2.23) 



According to Eq. ( |2.14| ) we then have a velocity field 

H z x r 



m r 2 



(2.24) 

in •/— 

which leads to a vortex energy 



2 



K=*pi— MR/Zo) + C (2.25) 



when inserted into Eq. ( 2.15 ), where, R is the system size. We have imposed a lower 
cutoff £o and included the contribution from smaller scales where the amplitude varies in 
the constant C. Since there are approximately (i?/£o) 2 distinct positions for a vortex, the 
free energy E v — TS V for large R is 

F v « (tt P °^ - 2Tj HR/Co) ■ (2.26) 

The prediction of this famous argument [|1]-[3| is that it becomes favorable for vortices to 
proliferate (leading to exponential decay of G(r)) above a critical temperature 

r e «~p2. (2.27) 

2 m z 



C. Vortex Statistical Mechanics and Renormalization of p s 

Going beyond the simple argument presented above requires a more detailed theory, 
which allows for a finite density of interacting vortices near T c . Understanding superfluidity 
in two-dimensional helium films is equivalent to solving the statistical mechanics of this 
vortex gas. We start by noticing that functional integrals like the partition function (|2.10 ) 
are dominated in the low-temperature "phase only" approximation by solutions of 

V 2 9 = . (2.28) 

To allow for vortices (and thus for amplitude fluctuations), we require that Eq. (2.28) be 



satisfied "almost everywhere," i.e., everywhere except in the cores of a collection N of vortices 
located at positions {r a } with integer charges {s a }. A vortex singularity has "charge" s a if 
the line integral of the phase gradient on any path enclosing the core satisfies, 

V0 • dt = 27ts a . (2.29) 
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The function ( |2.23| ) obeys this condition, for example, with s a = ±1. More generally we 
expect, for a contour C enclosing many vortices, that 

<£v6-di= f d 2 rn v (r) , (2.30) 
Jc Jn 

where Q is the area spanned by C and the vortex "charge density" is 

N 

n v (f) = 27r s a S(f- f a ) . (2.31) 

a=X 

Equation (|2~30| ) is a statement about the noncommutivity of derivatives of 9(x,y). In- 
deed, by taking the contour C to be a small square loop, we readily find 

e ij d i dj9(f) = d x d y 9 - d y d x 9 

= n v (f), (2.32) 

where is the antisymmetric unit tensor in two dimensions, e xy = —e yx = 1. To cast this 
equation in a more familiar form, we introduce the Cauchy conjugate to the phase field, by 

0i0(r) = eyfi^f) , (2.33) 

and find that 0(x,y) satisfies 

V 2 6(r) = n v (r) . (2.34) 

Finding the phase associated with a set of vortex singularities thus requires that we first 
determine the "electrostatic potential" 9(f) of a collection of point charges {sj} at positions 
{fj} in two dimensions IJ. The solution of Eq. ( |2.34j ) is 

9(f) =2irJ2 s a G(f,f a ) , (2.35) 

a 

where the Green's function satisfies 

V 2 G(f,f Q ) = 5(f-f a ) . (2.36) 
For \f — fj\ large and both points far from any boundaries, we have 

G(r, f 3 ) « i- In + C , (2.37) 

where C is a constant which contributes to the vortex core energy. 

We now decompose the phase into a contribution 9 v (x,y) from vortices, obtained by 
taking the Cauchy conjugate of Eq. Q2.35 ), and a smoothly varying part (fi(x,y), 
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9(x,y) = 8 v (x,y) + <f>{x,y) 



(2.38) 



The function <j)(x, y) represents single-valued phase fluctuations superimposed on the vortex 
extrema. Thermodynamic averages are obtained by first integrating over this nonsingular 
phase field, and then summing over all possible complexions of vortex charges and positions. 
To insert this decomposition into Eq. Q2.13| ), we need 



W(f) = 2n(z x V) J d 2 r'n(f')G(f,f') + V0 . (2.39) 
The resulting free energy takes the form |J 

£ = cons t. + l Ko J d 2 r\ V0| 2 + ^ , (2.40) 



where the vortex part is 



^ = -nK £ s a s, In (^-^) + ^ £ 4 > (2-41) 

1 a^f3 \ a / 1 a 

and the core energy E c is usually assumed to be proportional to K . Implicit in the statistical 
mechanics associated with Eq. fl2.41Q is a constraint of overall "charge neutrality," J2j &j — 0, 
required for a finite energy in the thermodynamic limit. 



D. Renormalization Group and Universal Jump in the Superfluid Density 

To illustrate the statistical mechanics of the vortex gas described above, consider the 
renormalized superfluid density pf(T) calculated to lowest order in the vortex fugacity 

y = e~ E ^ T . (2.42) 

The renormalized superfluid density is related to the correlations of the momentum density 
g(r). On a microscopic level, the momentum density operator gi{r) is given in terms particle 
creation and destruction operators by 

h 



= Ti 



ip + {r)diip(r) -^{r)diip + {r) . (2.43a) 



When local off-diagonal long-range order is present in helium films, we replace ip + (r) and 
ip(f) by coarse-grained two-dimensional classical c-number fields, as usual. In the "phase 
only" approximation which led to Eq. ( |2.13| ), the coarse-grained momentum density is then 



9{r) 



p%(r) , (2.43b) 
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where p° s = m\if)o\ 2 , and v s (f) includes possible contributions from vortices. In helium films, 
the contribution of the normal fluid to the momentum vanishes, due to the viscous coupling 
to the substrate. 

The correlation matrix which determines the renormalized superfluid density is 

C ij {q,K,y) = {g i {q)g*{q)), (2.44a) 

where gi{q) is the Fourier transform of the i-th component of g(f). To extract the renormal- 
ized superfluid density, we first decompose into transverse and longitudinal parts: 

Cij = A^f-f + B(q) Uj - . (2.44b) 

In an isotropic classical liquid, one would have A(q) = B(q) in the limit q — > 0. The mo- 
mentum fluctuations, moreover, would decouple from the configurational degrees of freedom 
responsible for most phase transitions. The behavior of quantum fluids is different: The 
renormalized superfluid density pf (T), in particular, is given by the difference between A 
and B as q tends to zero (see Appendix) 

pf (T) = ^ lire v [A(q) - B(q)} . (2.45) 
1 q^o 

The momentum density g(r) = p s —V6{f) is already decomposed into transverse and 



longitudinal parts in Eq. ( |2.39| ). One readily finds that Cy takes the form ( |2.44b[ ), with 



h 2 

m 2 T 



A(q) = K (2.46a) 



h 2 Att 2 K 2 

' -B(q) = —±(n v (q)n v (-q)} . (2.46b) 



m?T q 



where n v (q) is the Fourier transform of the vortex density n v (f), and the average in ( |2.46b|) 
is to be carried out over the vortex part of the free energy. 

When the vortices are dilute, it is straightforward to use these results to obtain a 
fugacity perturbation expansion for the renormalized superfluid density. In terms of 
K R = h 2 pf/m 2 T, we have § 

K^ = K^ + AnYo ~ ~) +O(y 4 o)- (2-47) 



a \a 



At low temperatures, Eq. ( |2.47| ) provides a small correction (proportional to e~ 2Ec ^ T ) to K^ 1 . 
When Kq ^ 2/ir, however, the integral becomes infrared divergent and perturbation theory 
breaks down. This is precisely the condition ( |2.27| ) required for a vortex unbinding transition. 
This potentially divergent perturbation theory can be converted into renormalization group 
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recursion relations for effective couplings K{£) and y(£) describing renormalized vortices 
with effective core diameter ae~ t . These differential equations read 

dK-\£) 



dy(£) 
di 



[2-7rK(£))y(£) + 0[y l 



(2.48a) 
(2.48b) 



Another important result is the invariance of the superfluid density along a renormalization 
group trajectory, 



K R (K,y) = K R {K{£),y(£)) 



(2.49) 



where K R (K, y) = -^pf(K, y). This result also follows from the Josephson scaling relation 
for pf E3. 




FIG. 4. Renormalization flows arising from the Kosterlitz recursion re- 
lations. The shaded domain of attraction of the fixed line at y(l) = is 
a superfluid. A locus of initial conditions is shown as a dashed line. The 
superfluid phase is bounded by the incoming separatrix which terminates at 
K~ l = 2/vr. 

Equations (2.48) are the famous Kosterlitz recursion relations, originally derived by him 
using another method Q. The Hamiltonian trajectories they generate in the (K , y)-plane 
are shown in Fig. 4, together with a temperature-dependent locus of initial conditions, 



-cK 



(2.50) 



where c is a constant. 
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Initial conditions to the left of the incoming separatrix renormalize into the line of fixed 
points at y — 0, which describes the low-temperature phase. At higher temperatures, y(£) 
eventually becomes large, indicating that vortices unbind at long wavelengths. One then 
expects that order parameter correlations decay exponentially, as in Eq. ( [2.21 ), where the 
correlation length £ + is related to the density of free vortices nf. 

n f (T)^^ 2 (T). (2.51) 

A variety of detailed predictions for the Kosterlitz-Thouless transition follow from these 
recursion relations and the transformation properties of various correlation functions under 
the renormalization group We focus here on the prediction of a universal jump discon- 
tinuity in the superfluid density [p3 |. This is a direct consequence of the relation (|2.49| ) and 



the renormalization group flows shown in Fig. 4. Below T c , y(£) tends to zero for large £, 
and we have 

K R (K,y)= lim K R (K(£),y(£)) 

= lim K(£) (T < T c ) . (2.52) 

Since it is clear from Fig. 4 that this limit is just 2/V, at T c , we have 

f) 2 o R (T) 2 

lim K R (K,y) = lim = - , (2.53) 

independent of the way in which the initial locus crosses the incoming separatrix. In helium 
films of varying thickness on varying substrates, one predicts a sequence of curves with jump 
discontinuities and T c values all falling on a line with slope 

pj£) = 2m 2 ^ 3 x iQ _ 9 ^ cm -2 K -i , 2M) 
T c tiTi k B T 
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slope = 3.491 x 10 gm/crn -°K 



.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 

T C (K) 

FIG. 5. Jump discontinuities in the superfluid density vs. temperature 
for over 70 different experiments for helium films. One representative p s (T) 
curve is shown. Adapted from Ref. [11|. 



Fig. 5 shows the locus of jump discontinuities in the superfluid density vs. temperature 
for over 70 experiments on helium films. All points lie on a line with slope close to the 
universal value (|2.54 ). This remarkable universal jump is related to the universality of the 
critical exponent rj(T) at the Kosterlitz-Thouless transition. When vortices are taken into 
account, G(r) 



i 



-^jj for T < T c , provided Eq. (|2.20|) is replaced by 
V (T) = l/2nK R (T) , 



(2.55) 



so that (neglecting logarithmic corrections) r)(T c ) = 1/4. 



E. Two-Dimensional Superconductors 

Superconducting films are similar to substrates coated with superfluid He 4 , with the 
understanding that the complex order parameter ip(r) describes Cooper pairs of electrons. 
Because the bosonic degrees of freedom are now charged, the currents surrounding vortices 
are screened by a coupling to the vector potential. The interactions between vortices are 
no longer logarithmic at all distances, but instead die off rapidly at scales larger than the 
London penetration depth. We first review a "naive" theory of screening in two-dimensional 
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superconductors. This theory is "two-dimensional" in the sense that all quantities are inde- 
pendent of the z coordinate. It would be directly applicable to situations in which vortices 
are, in effect, infinitely long rigid rods in a bulk material, and no thermal variation in the 
phase or vortex positions allowed along the z axis. The conclusion of the naive theory is 
that, due to screening, vortices will always unbind at any nonzero temperature fl[]. We 
then discuss real superconducting films, which behave differently, due to the spreading of 
the vortex magnetic field lines (see Fig. 6) as they emerge from the top and bottom of the 
film p4| . As emphasized by Beasley et al. |25| vortices in sufficiently thin superconducting 
films now interact logarithmically out to an effective London penetration depth of order 
millimeters, and the Kosterlitz-Thouless theory for neutral superfluids becomes directly ap- 
plicable over a wide range of length scales in most samples. The detailed predictions of the 
Kosterlitz-Thouless theory when transcribed to superconducting films are discussed in Ref. 



Supercurrents 



r 
B 




FIG. 6. Schematic of currents and magnetic field lines for a vortex in a 
superconducting film. Currents only flow in a two-dimensional plane, so the 
screening is weak. 



We begin with the "naive" theory of two-dimensional superconductors because it pro- 
vides a simple illustration of the screening of defect interactions by a vector potential. A 
closely related, but nonlinear screening phenomenon arises for dislocations in membranes 
(see Sec. Ill), and leads to the striking conclusion that crystalline membranes must always 
melt into hexatic liquids at any finite temperature. 
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1. Naive Theory 



We start with Ginzburg-Landau theory of superconductivity [EE], which generalizes Eq. 
( p. 11| ) for neutral superfluids to charged Cooper pairs, 



d 2 r 



2m* 



i c , 



1> 



+ _ r |^|2 +M |^|2 



87T 47T 



^2.56) 



Here, ip(r) is a coarse-grained complex field which represents the microscopic Cooper pair 
destruction operator (ip*(r) is the coarse-grained pair creation operator), A(r) is the vector 
potential, e* = 2e is the pair charge and m* = 2m is the effective mass of the pair, if is a 
constant external field, and B(r) = V x A(r) is fluctuating local magnetic field. Equation 
( p.56| ) is incorrect for real superconducting films because it implicitly assumes that both ip{r) 
and A(r) vanish outside the thin slab occupied by the superconductor. This is a reasonable 
assumption for ip(f), but incorrect for B{f) and A(r). This theory does, however, describe 
a bulk superconductor which is infinitely stiff in the z-direction, so that all quantities are 
independent of z. Equation fl2.56| ) should then be interpreted as a free energy per unit 
length. A related theory provides the correct description of local order in two-dimensional 
smectic liquid crystals |27[ . 

Let us see what Eq. ( |2.56| ) predicts for superconductors, with the external magnetic field 
H set to zero. The modifications introduced by the fluctuating magnetic field B(f) outside 
the sample will be discussed in the next subsection. Following our treatment of helium 
films we assume that thermal excitations have suppressed any possible ordering well below 
the mean field transition temperature, so that r C 0. We again make the "phase only" 
approximation, and write 



(2.57) 



with ipo = ^J—r/Au. The two terms quadratic in ip(r) determine the superconducting 
coherence length, 



£o = \J h 2 / m*\r\ 



(2.58) 



which is the length over which the order parameter rises from zero near a vortex. On scales 
larger than £ 0) so that the phase only approximation is valid, free energy becomes 
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where 



^ = m*|Vor (2-60) 

is the mass density of the Cooper pairs. 

To see qualitatively how the vector potential affects the behavior of vortices, imagine that 
the vortex solution ( |2.23| ) for the phase variable is inserted in Eq. ( |2.59|) . Since |V6*| ~ 1/r, 
this would lead to a logarithmically diverging energy, if the vector potential were neglected. 
The vector potential, however, will cancel out this divergence and reduce the energy provided 
it too falls off like 1/r with the correct coefficient far from the vortex. This screening costs 
gradient energy, however, due to the |V x A\ 2 term. The two terms of Eq. ( |2.59| ) quadratic 
in A become comparable when A(r) varies over a length scale Ao, with 

(2.61) 



Ao 



(m*c) 2 



Ti(e*) 2 



This scale is the London penetration depth, of order 1-10 thousand Angstroms in most 
superconductors. Because of this is the scale over which screening sets in, we expect the 
vortex energy to be finite and approximately equal to Eq. ( |2.25| ) with R replaced by A , 

E v n*p° a (J^ ln(cA /£o) (2.62) 

where c is an undetermined constant. Since the entropy of an isolated vortex is still pro- 
portional to lni?, and the energy is now finite, we conclude that vortices will always be 
unbound in this model [|I|]. A more detailed analysis, similar to that for finite energy defects 
like vacancies and interstitials in a solid ||28|| , predicts a nonzero density of free vortices 

Vf « Ho 2 e~ Ev/T (2-63) 

at any temperature T > 0. 

Because Eq. Q2.59Q is a simple quadratic form, we can go beyond the arguments sketched 



above, and determine the interactions between screened vortices in more detail. As in helium 
films, functional integrals weighted by exp[— F/T] will be dominated by vortex extrema of 
the free energy. In the Coulomb gauge, V • A = 0, the variational equation 8F/59 = again 
leads to V 2 6 = away from the vortex cores, and we can immediately write down the vortex 
contribution to V6 from iV vortices at position {f}} with "charges" {sj}, 
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W(f) = 2vr s a (z x V)G(f, r ) , (2.64) 



where G(f, r}) is the Green's function discussed in Section II. C above. Note from Eqs. (|2.31 
and ( |2.32| ) that the curl of V6 1 does not vanish, 

z- (V x V0) =2nJ2sJ {2) (r-r a ) . (2.65) 

a 

The vector potential is constrained by vortex singularities, even though 9{r) and A(r) 
appear to decouple in the Coulomb gauge. To see this, consider a closed path C bounding 
an area Q in the superconductor far compared to Ao from any vortex cores. On such a path 
the first term in ( 2.59Q has been minimized, insuring that V6 1 = and that 



. (ve-^rA) -d£ = 0. (2.66) 

c \ he ' 



It then follows from Eq. ( |2.3U| ) that the magnetic flux through VL is given by the sum of the 
enclosed vortex "charges," 



d 2 rB z = <j) ]T s a , (2.67) 



where 0o — ^fr^ is the flux quantum. 

A second variational equation follows from 5F/5A = 0, i.e., 



V x [V x A(f)] = V x B(f) = Anp° s \JL\ (v6 - Q A(r) . (2.68) 

Upon casting this equation in the form of Ampere's law V x B(r) = 4irJ(r)/c, we see that 
the (gauge invariant) supercurrent associated with 9 and A is 

J(f) = e^ | 2 A (W(f) - £ A(r)) . (2.69) 
m* V he ) 

Note the close analogy with the formula for the superfluid momentum density in helium 



films 



g{f) = m\ipo\ 2 —V8(f) . (2.70) 
m 

A closed form equation for the magnetic field B z generated by a distribution of vortices 
follows from taking the curl of Eq. (|2.68|) , and using Eq. (|2.65| ) 

N 

B z - \lV\B z = 0o E sj {2 \r- r a ) , (2.71) 

a=l 
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where "_L" denotes coordinates in the plane of the film. 

Consider one isolated vortex with charge s a at the origin. The Fourier transformed field 
B x (q±) which solves Eq. fl2.71|) is then 



S(p 



BM = j^j , (2.72) 



which leads in real space to 

S<f)Q 



B *(r) = ^T2^o(r/A ) , (2.73) 



where K (x) is the Bessel function, K (x) ~ ln(l/x), x -C 1, and K (x) = ^Jtc/2x e x , 
x ^> 1. The associated supercurrent J(r) = V x B(r) is 

OVT Ag 

where is a unit vector in the azimuthal direction. Note that \J{r)\ ~ 1/r for r Ao, 
similar to the momentum density in helium films. However, screening sets in for r ^> X Q , 
and | J(r)\ ~ exp(— r/Ao). Similar manipulations lead straightforwardly to the generalization 
of the vortex-free energy (|2.41 ) for two-dimensional charged superfluids, 



i^ = eoE s ^i^Q + ^ E s ? > ( 2 - 75 ) 



where eo = (0o/4vrAo) 2 and S c is a vortex core energy. Gaussian "spin wave" fluctuations 
about these vortex extrema can also be included [|27|. The finite range of the interaction 
potential confirms our earlier conclusion that vortices will unbind for entropic reasons at 
any finite temperature in this model. 



2. Real Superconducting Films 

As illustrated in Fig. 6, the magnetic field lines generated by a vortex in real super- 
conducting films spread out into the empty space above and below the plane of film. This 
spreading leads to different behavior for the screening currents in the film than that predicted 
by Eq. (|2.74j ). The true behavior was first found by Pearl ||24j| , and provides an interesting 



illustration of how two-dimensional physics can be effected by the outside, three-dimensional 
environment. We sketch the calculation here, following the treatment of deGennes p0 |. 



Start by integrating Eq. ( |2.69D over a film of thickness d along z. Provided the film 



is much thinner than the London penetration depth, we may assume that 8(f) and v4 2 d(r) 
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remain constant across the film thickness. We denote the vector potential by A 2d (r) to 
emphasize that this quantity is the vector potential in the two-dimensional film, and not the 
full vector potential A 3d (f, z) which describes the magnetic field in all of three-dimensional 
space. Note that r = (x, y) always refers to coordinates in the two-dimensional plane 
perpendicular to z. If we imagine that all this current is concentrated in a delta function 
sheet in the plane z = 0, the full three-dimensional current may be written 



/(?,*) c ^ 



An 27rA e ff 
where 



2ir -. 
W(r) - —A 2d (r) 

00 



6{z) , (2.76) 



A eff = \ 2 /d . (2.77) 

As we shall see, A e fr will play the role of an effective London penetration depth in this 
problem. 

The current (|2.76|) provides a source for a i/iree-dimensional magnetic field B(r,z) via 
Ampere's law, V x B(r, z) = 4irJ(f, z)/c. Upon setting B = V x and using the Coulomb 
gauge, Ampere's law becomes 



- V 2 A 3d (r, z) + —A 2d (r)5(z) = —5(z) , (2.78) 

where V 2 = + d 2 /dz 2 and we have inserted the phase gradient for a single vortex of 
charge s at the origin. We now pass to Fourier transformed vector potentials 

&d(q±, Qz) =J d2r J dze^e^ z A 3d (f, z) (2.79) 

and 

A 2d (q ± ) = J d 2 re l ^A 2d (f) (2.80) 

and find that (12.781) becomes 



7 t- 1 7 /- is<Po z x Q± 

A 3d {q±, q z ) + , 2 2 , A 2d (q ± ) = — . (2.81) 

Ks(q 2 + q±) KM + q±) q± 

The 2d and 3d vector potentials must agree on the plane z — 0, i.e., A 3d (r,z = 0) = 
A 2d (f), so that 

A 2d (q±)= f°° ^A 3d (q±,q z ) ■ (2.82) 
After integrating Eq. (|2.81|) over q z and using this result we can solve for A 2d (qj_), 
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A 



2d{q±) 



ts(p {z x q ± ) 
ql(l + 2X cS q ± ) 



(2.83) 



Upon inserting Eq. Q2.83 ) into the Fourier transform of Eq. ( |2.76| ) after setting J(r, z) 
dJ2d(r)S(z), we have 



J2d{qi 



c is<po{z x q ± ) 
2n dqj_(l + 2X eS q ± 

and so the superconducting current in real space is 



2nd 



-iqj_-r 



z x q± 



(2.84) 



(2.85) 



(2tt) 2 (1 + 2X eS q ± ) q ± 

The behavior of this screening current depends on the effective London penetration 
depth. For r <C A c g, we find 

SC0Q 



J 2 d(r) 



8n 2 dX e f[r 







(2.86) 



which agrees with the r Ao limit of Eq. (|2.74|) for our "naive" model. This unscreened 
1/r falloff now continues, however, out to a much larger distance A e fj = (A /<i)Ao ^ A . For 
t A c flf, the limiting behavior of (|2.85| ) is 



SC(p 







(2.87) 



An 2 dr 2 

The exponential screening of vortex currents in the "naive" model is thus replaced by a 1/r 2 
power-law falloff. 

Similar results hold for the interactions between vortex pairs [0. Vortices interact 
logarithmically, as into neutral superfluids for r -C A e fr, but exhibit a weaker 1/r potential 
for r ^> A fj. Since A e g can be of order fractions of a centimeter for film thicknesses d = 
10 — 100 A, the Kosterlitz-Thouless theory becomes directly applicable on essentially all 
length scales for sufficiently thin films 



. The relatively weak screening compared to bulk 
systems arises because currents are confined to a thin plane instead of forming rings along 
the entire z axis. 



III. DEFECTS IN MEMBRANES AND MONOLAYERS 

Two-dimensional crystals have much in common with two-dimensional superfluids. There 
are many important experimental examples, including rare gases adsorbed onto periodic 
substrates like graphite, Langmuir-Blodgett films of ampiphillic molecules at air-water in- 
terfaces, freely suspended liquid crystal films, electron layers trapped at the surface of liquid 
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helium, and assemblies of colloidal particles confined between two glass plates P,|31fl. Crys- 
tals consisting of a few atomic or molecular layers display algebraic decay of a translational 
order parameter, similar to Eq. ( [2.19| ), and the crystal elastic constants play a role similar 
to that of the superfluid density || . As emphasized by Kosterlitz and Thouless [|IJ] and by 
Berezinski [0, dislocations in such crystals are point defects with a logarithmically diverging 
energy as a function of system size, and one might expect them to melt via a dislocation 
unbinding mechanism. 

When the detailed dislocation unbinding theory was worked out |]32|j33| , there was a 
surprise. Melting via dislocations leads not to an isotropic liquid, as proposed originally 
mill, but produces instead a new hexatic phase of matter, with residual bond orientational 
order [32]|. The long-range bond orientational order in two-dimensional crystals at low 
temperatures is converted into algebraically decaying correlations in a hexatic fluid by a 
gas of unbound dislocations. Each dislocation, moreover, contains an embryonic pair of 
orientational defects called disclinations in its core. These disclinations separate and interact 
logarithmically in the hexatic phase, and themselves unbind via a second phase transition at 
sufficiently high temperatures. The latent heat associated with the usual first-order melting 
point can thus be spread out over an entire intermediate phase, separated from the low- 
temperature crystal and high-temperature liquid by two continuous phase transitions. 

We call the experimental systems mentioned above "monolayers," to emphasize that they 
are constrained to be approximately flat. In all cases, there are nearby substrates, walls or 
interfaces which force the degrees of freedom to "layer," i.e., to lie a plane. Excitations out 
of the plane are strongly disfavored by, for example, a surface tension. There is, however, an- 



other important class of two-dimensional materials called "membranes" [34]. Membranes are 
two-dimensional associations of molecules different from the three-dimensional fluid medium 



in which they are embedded. Examples include lipid bilayers in water [33] or spectrin pro- 
tein skeletons extracted from red blood cells ||35|| . Because membranes are not confined to 
an interface between two different phases, the surface tension vanishes and they exhibit wild 
fluctuations out the plane while retaining a local two-dimensional topology. 

In this section, we review the physics of defects such as dislocations, disclinations, grain 
boundaries, vacancies and interstitials in membranes and monolayers. We discuss crystalline 
monolayers first, and then show how defect buckling screens defect energies in membranes 
PB^pS] . The buckling transition of disclinations in hexatic membranes is also described. 
We point out a fundamental asymmetry in the populations of positive and negative discli- 
nations in liquid membranes with free boundary conditions or with fluctuating topologies. 
For a discussion of the spin-glass-like statistical mechanics of polymerized membranes with 
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quenched random distributions of defects such as impurity atoms, grain boundaries, see 
Refs. HHH. 



A. Landau Theory and Elasticity of Tethered Membranes 

Consider first the statistical mechanics of two-dimensional assemblies of atoms and 
molecules without defects. To exclude defects, we assume a perfect triangular lattice of 
monomers embedded in three dimensions and tethered together via an unbreakable network 
of covalent bonds. Such "tethered surfaces" [4^,44] exhibit interesting fluctuations and phase 



transitions even in the absence of defects. Once in the low temperature, broken symmetry 
phase of the relevant Landau theory, defects can be introduced for both monolayers and 
membranes. 

Consider first the high temperature, crumpled phase of a tethered surface, similar to 
the crumpled state of a linear polymer chain. We describe the membrane by a function 
f(xi,X2), where f is a three-dimensional vector specifying the position of the monomers 
as a function of two internal coordinates x\ and X2 fixed to the monomers. It turns out 
that the crumpled membrane can undergo a spontaneously symmetry-breaking into a flat 
phase below a crumpling temperature T c []3q| . The order parameters for this transition are 
the surface tangents {t a = dr/dx a , a = 1,2} and the Landau free energy describing the 
transition is ||45|| , 

1 



F[r\xi, x 2 )] = J d x 



+ -^o-idar) 2 + b(d a f- dpr) 2 + c(dyf ■ dyff + 



(3-1) 



Within this expansion in the tangents, the probability of a surface configuration f(xi,X2) 
is proportional to exp{— F[f(xi, X2)]/T}. Note that F is invariant under translations and 
rotations both within the embedding space and in the internal coordinates, as it should 
be. Self-avoidance between monomers at very different values of x = (xi,^), but close 
together in the embedding space is neglected, although it is not hard to incorporate this 
into the model [ |45|| . We assume that a changes sign at the mean field crumpling transition, 
a = a'(T — T c ), just as in the more conventional Landau theories for superfluids discussed 
in Sec. II. 

We discuss here only the low-temperature flat phase. When a < 0, F is minimized by 

f = (f(x 1 ,x 2 )) 

= m\x x ex + x 2 e 2 } , (3.2) 
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where e\ and e 2 are an arbitrary orthogonal pair of unit vectors in the three-dimensional 
embedding space and 



m 



2 V b + 2c 



(3.3) 



Note that (t a ) = (df/dx a ) = me a , so that m is the "amplitude" of the tangent order 
parameters. The average of the surface tangents vanishes in the high-temperature phase, 
because the surface is highly crumpled. Additional physical insight is provided if we rewrite 
the free energy as 

1 2 



F = const. + J d 2 x |^((9^f) 2 + /i 

+ A 

where 
and 
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(3.4) 



(3.5a) 



(3.5b) 



This free energy is a sum of a bending energy, controlled by k, and stretching contributions 
governed by elastic moduli fi and A. The stretching terms provide an energetic penalty 
whenever the induced metric in the embedding space, 

df df 



g a p{xx,x 2 ) 
deviates from a flat background value 



dx Q dxp 



g°p = m 2 5 aP . 



(3.6) 



(3.7) 



The low energy Goldstone modes associated with the flat phase are phonons. To study 
these excitations, we proceed as in the "phase only" approximation for superfluids j|(J and 

set 



f(x , x ) w m(xi + «i)ei + m(x 2 + u 2 )e 2 + fe 3 



(3.8) 



where u(xi,x 2 ) is an in-plane phonon field and f(xi,x 2 ) is an out-of-plane displacement 
along e*3 = vece± x e 2 . With the neglect of an additive constant, the free energy 
becomes 



27 



F = J d 2 x\^[d a up + d(,u a + {d a f){dpf)} 2 

+ \[d,u 1 + \{d 1 fff + 'd 2 x{V 2 ff • (3.9) 

This expression is identical to the energy of a bent elastic plate, where k is bending rigidity, 
fi is a shear modulus, and fi + A is the bulk modulus pETfl . One often finds the stretching 
energy contributions written in terms of the nonlinear strain matrix w«,«(£i, X2) 

Ua(3 = ^(d a up + d p u a ) + ]^{d a f){dpf) . (3.10) 

We have neglected terms nonlinear in u, which are less important than the nonlinear terms 
in / which we have kept. It can be shown that the nonlinear terms in / cause the renormal- 
ized bending rigidity to diverge at long wavelengths, leading to the remarkable conclusion 
that the low-temperature flat phase represents a genuine broken continuous symmetry with 
true long-range order [pi|p6[ . Broken continuous symmetries two-dimensional systems at 
finite temperatures are usually impossible, and are replaced instead by exponential decay 
of correlations or by algebraically decaying order such as that found for superfluid helium 
films in Eq. ( g!9|) . 

Note the similarity between Eq. (|3.9|) the free energy Eq. (|2.59|) which arose in our 
"naive" theory of superconducting films. The phase gradient V# is replaced by the linearized 
strain matrix \(d a up + dpu a ). The superfluid density is replaced by the elastic constants 
fi and A. The role of the "vector potential" is played by gradients of /, and the bending 
rigidity contribution replaces the field energy term g^|V x A\ 2 for superconductors. Because 
"vector potential" contribution to Eq. ( |3.10| ) is nonlinear, the two theories are certainly 
not identical. Nevertheless, we shall find striking similarities when we relax the constraint 
of perfect sixfold coordination and allow for defects: If we first set / = 0, then Eq. ( |3.9| ) 
can be used to calculate defect energies in planar "monolayers" — see below. Defects in 
monolayers are like the vortex singularities discussed earlier for neutral superfluids. When 
/ 0, however, the theory applies to the low-temperature phase of "membranes." As we 
shall see, these defects are then screened by nonzero gradients of /, similar to the screening of 
vortices by the vector potential in superconductors. The corresponding reduction in energy 
for membranes is accompanied by defect buckling out of the plane defined by / = 0. 
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B. Defects in Monolayers 



Fig. 7 shows a typical particle configuration for a thermally excited monolayer in its crys- 
talline phase. Most particles, indicated by circles, have six nearest neighbors, as determined 
by the Dirichlet or Wigner-Seitz construction: The nearest neighbor coordination number of 
a particle is the number of edges of the minimal polygon formed by the bisectors of the lines 
connecting it to its near neighbors. The 5- and 7-fold coordinated particles (indicated by 
diamonds and asterisks, respectively) are orientational disclination defects in the otherwise 
6-fold coordinated triangular lattice. An elementary application of Euler's theorem shows 
that the average coordination number with periodic boundary conditions must be exactly 
six [48|, so the numbers of 5's and 7's must be equal. A high-temperature liquid can be 



viewed as a dense plasma of disclinations. The "plus" and "minus" disclinations (5's and 
7's) in a liquid annihilate and pair up with decreasing temperature to form a crystal. 
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FIG. 7. Snapshot of a particle configurations in a computer simulation 
of a "monolayer." Six-fold coordinated particles are shown as circles. Five- 
fold and seven- fold disclination defects are shown as diamonds and asterisks, 
respectively. Heavy contours surround a dislocation and an interstitial. From 
a simulation of classical electrons interacting with a repulsive 1/r potential 
by Rudolph Morf. 
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Disclination pairing is very evident in Fig. 7. An isolated disclination pair is, in fact, a 
dislocation defect. As illustrated in the figure, dislocations are characterized by the amount 
a path which would close on a perfect lattice fails to close. This mismatch, or "Burgers 
vector," is a lattice vector of the underlying triangular crystal. It acts like a discrete vector 
"charge" attached to the dislocation, is independent of the exact contour chosen and points 
at right angles to a line connecting the 5 to the 7 in the dislocation core. More generally, 
such circuits determine the vector sum of all dislocation Burgers vectors contained inside. 
The Burgers circuit does in fact close for the other circuit shown in Fig. 7. The necklace of 
alternating 5's and 7's it contains can be viewed as three dislocations with radial Burgers 
pointing at 120 degree angles to each other. It represents a third type of defect, an interstitial. 
To see this, note that since the closed Burgers parallelogram is 4 lattice units on a side, it 
would enclose 3x3 = 9 particles if the lattice were perfect. In fact, there are 10 particles 
inside the contour, indicating the presence of an extra atom, or interstitial. The presence of 
a vacancy, the "antidefect" of the interstitial, could be detected by similar means. 

To calculate the energies of the various monolayer defects discussed above, we use the 



continuum elastic free energy (3.9) with / = 



F=\J d 2 r[2^u% + \u 2 kk ] (3.11) 

where 

u ij = + djUi) . (3-12) 

In experimental monolayer systems, the constraint / = might be imposed by a strong wall 
or substrate potential. For particles confined at an interface with a surface tension a, we 
should add a term 

SF = la J d 2 x{Vf) 2 (3.13) 



to Eq. (313). Upon integrating out the / field in perturbation theory, and we recover a free 



energy of the form ( |3.11| ), with renormalized elastic constants ji and A ||49|| . 



Dislocations, disclinations and other defects can be introduced into the theory in a way 
similar to the discussion of superfluid vortices in Sec. II. First, however, we determine the 
equations satisfied by free energy extrema away from defect cores. The variation of the 
monolayer free energy ( |3.11| ) with respect to u leads to 

0^=0, (3.14) 
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where 



dij = 2/iUij + XukkSij (3.15) 

is the stress tensor. Equation ( |3.14j ) will be satisfied automatically if we introduce the Airy 
stress function \ Yia - 

Vijir) = £imejnd m dnX(r) . (3.16) 
The individual components of are thus 

d 2 X d 2 X d 2 x m? x 

° xx dy 2 ' ayv dx 2 ' Gxy dxdy ' 1 } 

The function x(r) is similar to the vector potential one uses to insure V • B = in Maxwell's 
equations. If we are able to find Xi we know cr^, and hence the strain matrix by inverting 
Eq. QSTp , 



= 77- e imejn5 m 9 n x - "j — 7~~7~TT^ 2 X • (3.18) 
2/x 4/i(/i + A) 

So far, we have not used the important fact that Uij is determined by gradients of a 
displacement field via Eq. ( |3.12|) . The corresponding requirement on the superfluid velocity, 



as given by Eq. (|2.14j ), is 

Vx^O, (3.19) 

which must hold away from vortex cores. The analogous compatibility condition on 
follows from applying the operator e^e^c?^ to both sides of Eq. ( |3.18| ). This operator van- 
ishes when acting on iiij(r), provided various derivatives of the displacement field commute. 
Equation (|3.18|) then simplifies to give a biharmonic equation of x(^)> 



J_V 4 x(r) = ~eij£jid k di(diUj + djUi) 
K Z 

= , (3.20a) 

with 

0.20b, 
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(a) 




(b) 

FIG. 8. Five- (a) and seven-fold (b) disclinations in monolayers (flat) and 
membranes (buckled). From Ref. [37]. 

Equation (3.20) must be satisfied "almost everywhere," that is away from defect cores. 
Defects introduce source terms on the right-hand side, similar to the vortex density which 
appears in Eq. (|2.34 ). As for superfluid vortices, defects represent points where derivatives 
fail to commute. Because the line integral of the differential bond angle field 8(f) around a 
disclination in a triangular lattice must be an integral multiple of 2tt/6 = 60°, 



/ d8(f) = s — 
Jc 6 



(3.21) 



we have 



7T 



eijdidj9(r) = '-^s a 5(r- f a 



(3.22) 



for a collection of disclinations with "charges" {s a = ±1, ±2, . . .} at positions {fa}- Isolated 
5- and 7-fold disclinations are shown in Fig. 8. The local bond angle 8(f) in a crystalline 
solid is given by the antisymmetric part of the strain tensor, 

1 



8(f) = [d x Uy(f) - dyU X (f)} 



(3.23) 



so the noncommunitivity of derivatives for disclinations takes the form, 
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1 7T 

-eijtkedidjdkUi = - ^ s a 6{f - f 



(3.24) 




(a) 




(b) 

FIG. 9. Dislocation defect in a (a) monolayer (flat) and a (b) membrane 
(buckled). The burgers construction is shown in both cases. From Ref. [37]. 

We can treat dislocations by regarding them as tightly bound disclinations pairs, or 
alternatively by starting with the definition of the Burgers vector in the continuum limit 

dm 



r . , d,Xj 

C OXj 



(3.25) 



where (see Fig. 9) the circuit C encloses a single dislocation with Burgers vector b. For 
a collection of Burgers "charges" {b a } at positions {r a }, the differential statement of the 
noncommutivity embodied in constraints like ( |3.25| ) is 



(3.26) 



Upon noting that the right-hand side of Eq. ( |3.20a| ) may be rewritten, 

2 e ikej£d k de(diUj + djUi) = e M d k di>d + e ip e M d p did k Ui , 
we have, combining the above results for disclinations and dislocations, 

1 7T 

_V 4 x(r) = eijb a idj5(r - r a ) + ^Y,sp6(f- fp) . 

A a Op 



(3.27) 



(3.28) 
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Vacancies and interstitials are special cases of a general class of "impurity defects" which 
includes as well substitutional atoms of the wrong size. If the defect sits in a lattice site 
of 3-fold or higher symmetry, its stress field is isotropic, and it can be characterized by the 
area deficit Qq (positive or negative), it induces in an otherwise perfect lattice, 

n = /(V • u) d 2 r, 

h ■ u(r) di (3.29) 



ic 

where n is an outward unit normal to a contour C surrounding an impurity at the origin. 



The biharmonic equation for x( r ) reads [)47|1 



V 4 x(f) = -2/A)V 2 5(f) . (3.30) 

A similar equation can be derived for vacancies and interstitials by regarding them as three 
dislocations at the vertices of a small equilateral triangle with radial Burgers vectors pointing 
inward (vacancy) or outward (interstitial). 



It is straightforward to rewrite the free energy ( |3.11|) in terms of x( r )i 

1 



2K, 



d 2 r[V 2 X (r)] 2 . (3.31) 



Once x{^) is known, we can thus calculate both the strain fields (from ( j3.18|) and the defect 
energy. For an isolated disclination at the origin with angle defect s, we have pU 



and 



x(r) = —^—r 2 [\n(r /a) + const.] (3.32) 

07T 



K a 2 

03 R 2 , (3.33) 



327T 

for a circular patch of crystal with radius R. Because disclination strains do not fall off 
at large distances, the free energy diverges quadratically with the size of the system. We 
have neglected a core energy contribution which is negligible as R — ► oo. For an isolated 
dislocation at the origin, one has 

K 

X(r) = -r b i£ijrj\n(r/a) (3.34) 

and 

K b 2 

F = \n(R/a) + const. (3.35) 

87T 
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The logarithmic divergence arises because the strains fall off like 1/r, similar to the superffuid 
velocity in helium films. 

The energy of isolated disclinations is prohibitively large in crystalline monolayers. This 
energy is reduced to a logarithmic divergence in the hexatic monolayers, however, because 



of screening by a gas of unbound dislocations |32|] . The logarithmic energy of dislocations 
in two-dimensional crystals leads via the usual energy /entropy argument to an estimate of 
the dislocation unbinding temperature [0,0 . For a more detailed discussion of the statistical 
mechanics of defect mediated melting of monolayers, see Ref. ||. 

For completeness, we note that the Airy stress function of an isolated impurity at the 
origin is 

x (f) = Z^Ei n ( r / a ) (3.36) 

7T 

corresponding to a displacement field 

Because the strains now fall off as 1/r 2 , the elastic energy is finite and comparable to the 
core energy. 



u(r) = ^_ . (3.37) 



C. Defects in Crystalline Membranes 

We now return to the full continuum elastic free energy Eq. ( |3.9|) and determine how 
the defect energies are reduced in membranes, due to buckling into the third dimension. 



Variation of F with respect to / and u leads to the von Karman equations [47 



kV 4 / = + - 2-^2L-^L (3 38a) 

dy 2 dx 2 dx 2 dy 2 dxdy dxdy 

^ 1 ^-Uw + & +siF> - (3 ' 38b) 

The Airy stress function is related to the stress as for monolayers, (Jij(r) = €i m ej n d m d n x(r) , 
where we again have = 2fiUij + Xiikk&ij- Now, however, Uij is the nonlinear stress tensor, 

(3 - 39) 
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In the first line of ( |3.39| ), we have included nonlinearities in the in-plane displacements 
u. Although the nonlinearity in / is always important for membranes, the |J^- ■ term 
can often be neglected. This terms does contribute significantly to disclination energies, 
however: In monolayers, for example, it leads to a small correction to the coefficient of the 
R 2 divergence in the disclination energy for monolayers f50fl . In contrast to disclinations, 
the strain fields for dislocations and impurities fall off fast enough to justify neglecting the 
liar ' Jx 7 term in evaluating the far- field elastic energy. We have included a source term on 
the right-hand side of Eq. ( |3.38b|) due to defects. For an isolated defect at the origin, one 
has 

|s5(r), disclination, 
S{r) = I 6ijbidj5(f), dislocation, (3.40) 

> -^TX) fi oV 2 5(f), impurity, 

just as in flat space. The additional contribution to the right-hand side of Eq. (|3.38b|) 
is the Gaussian curvature associated with the out-of-plane membrane displacements. This 
"curvature charge" can partially cancel the "topological charges" S(f) of the defects when 
they buckle. Note that any solution with / 7^ 0, implies another solution with / — ► — /. 
When / = 0, we recover the elastic equations for monolayers. 



1. Disclinations and Dislocations 

Equations (3.38) are nonlinear and extremely difficult to solve analytically. Some 
progress is possible, however, in the inextensional limit Kq — > 00. All defect energies would 
be infinite in this case for monolayers. For membranes, however, defects can buckle so as 
to eliminate the elastic contributions to Eq. ( p.9|) . The only remaining contribution is the 
bending energy. The required / field is determined by setting the right-hand side of Eq. 
( p.38b| ) to zero. The inextensional solutions of the von Karman equations representing a 
5-fold disclination are 

x {r) = -Kln(r/a ), f(r) = ±-i/| r , (3.41) 

with energy 

E 5 = -nKhi{R/a ) , (3.42) 

for a circular membrane with radius R. A 7-fold inextensional disclination can be represented 
in polar coordinates approximately by [j36j , [37ll 
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2 

X(r) = 3Kln(r/a ), /(r) = ±\ - r sin 20 (3.43) 
which leads to an energy 

E 7 = TTKln(R/a) . (3.44) 

The R 2 divergence of disclination energies in monolayers (with an infinite coefficient when 
K — > oo!) is thus screened considerably by buckling. 

Precise numerical solutions of the von Karman equations for arbitrary k and Kq can be 



constructed by minimizing the energy of a triangulated tethered surface model p7 |. The 
bending rigidity is represented by an interaction between neighboring unit normals to the 
triangular plaquettes and the in-plane elasticity by nearest neighbor harmonic springs ad- 
justed to give the correct value of K . This "dynamic triangulation" method is easier to 
implement than a direct numerical solution of the von Karman equations, and includes 
automatically all relevant nonlinearities. The numerical results for disclinations in the inex- 
tensional limit are similar to Eqs. ( 3.42 ) and Q3.44 ), with slightly different coefficients [37 



It is still true that F§ < F-j. Buckled disclinations obtained using this approach are shown in 
Fig. 8. Five-fold disclinations buckle into a cone, while the seven-fold disclination leads to 
a saddle surface. For disclinations and dislocations, the crucial physics lies in the existence 



of a buckling radius |56 |: For small system radii R, these defects lie flat with no bending. 
Above a certain critical radius R c , the defects trade elastic energy for bending energy and 
buckle out-of-the-plane. Numerical studies [^] show that disclination energies are screened 
down to a logarithmic divergence beyond the buckling radius, as suggested by results in the 
inextensional limit. Both disclinations and dislocations behave as if they were inextensional 
beyond the buckling radius. 

Buckling is even more important for dislocations. Beyond the buckling radius Rb, the dis- 
location energy no longer increases logarithmically with R, and in fact appears to approach 



a finite constant | 37| . A buckled dislocation is shown in Fig. 9. The entropic contribution 
to the dislocation free energy = — SjT, however, still varies logarithmically with 
system size. The finiteness of the dislocation energy (or more precisely, any R dependence 
of the energy increasing more slowly than ln(R/a)) then implies that untethered crystalline 
membranes must melt at all nonzero temperatures [|36|j37f| . The areal density of dislocations 



in a hexatic membrane at temperature T is approximately 

n d « a^ 2 e~ Ed/T (3.45) 
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where Ed is is the (finite) dislocation energy, and ao is the lattice constant. Dislocations 
in membranes thus behave similarly to vortices in the "naive" model of superconducting 
films. Since buckled disclinations still have a logarithmically diverging energy, the hexatic 
fluid should be separated from a high-temperature isotropic liquid by a finite-temperature 
disclination unbinding transition (see below). Melting of membranes confined between two 
flat plates, which suppresses the buckling of dislocations and leads to a finite melting tem- 
perature T m (d), tending to zero as d — > oo, has been studied by Morse and Lubensky |H| . 



Defects like dislocations and disclinations buckle into the third dimension when the cost 
in stretching energy to remain flat exceeds the bending energy required to buckle. More 
generally, buckling occurs whenever 

^>7 (3.46) 



where i is a characteristic length scale for the defect and 7 is a dimensionless constant, 
typically of order 10 2 . In circular membranes of radius R, £ — \^Rb for a dislocation with 
Burgers vector b, and 7 ~ 127. For positive disclinations, we have I = R and 7 ~ 160, 



while for negative disclinations i m R and 7 ~ 192 [37]. Thus isolated dislocations and 



disclinations always buckle in sufficiently large crystalline membranes, irrespective of the 
values of the elastic constants. 

Strictly speaking, these conclusions apply only at T = 0. In a finite temperature defect- 
free flat phase, thermal fluctuations cause the long wavelength wavevector dependent bend- 
ing rigidity to diverge, K R (q) ~ 1~ Vk , while the renormalized elastic parameter K R (q) tends to 
zero K R (q) ~ q Vu |36|j52]j53|] . It is then appropriate to substitute K /k — > K r /k r ~ l/RP K+r,u 
in Eq. ( |3.46| ). Using the values r\ K = 0.82 and i] u = 0.36 p5| , we conclude that disclinations 
{£ ~ R) will still buckle, while dislocations (£ ~ R 1 ^ 2 ) become asymptotically flat as R — > 00. 
Dislocations could still buckle on short scales, before these long wavelength thermal renor- 
malizations of the elastic parameters set in. The dislocation energy will in any case tend to 
a finite constant for large R: Even if the dislocation remains asymptotically flat, we should 
now replace K by K R (R) ~ l/R Vu in Eq. (|3.35| ). Thus the conclusion that the hexatic 
phase is the stable low-temperature phase of membranes is unchanged. 



2. Other Defects in Crystalline Membranes 

Buckling is not inevitable even at T = for finite energy defects such as vacancies, 
interstitials, or impurity atoms. The buckling criterion is again Eq. (|3.46 ), where i is 
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now related to the excess area induced by the defects, I ~ yHo | 58| . Buckling can be 
triggered in an infinite system simply by varying the ratio K /k>. Because of the finite 
length scale associated with these defects, it is the bare "local" values of K and k which 
should appear in ( |3.46|) . Since k usually increases with increasing temperature, and K 
usually decreases, buckling is most likely at low temperatures. The / — > — / symmetry of 
the von Karman equations insures that buckled defects have at least two degenerate minima, 
representing displacements on opposite sides of the membrane. The ensemble of two-level 
systems generated by buckled finite energy defects will contribute to the specific heat and 
other equilibrium properties of membranes. 
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FIG. 10. Crushed vacancy defect in a flat monolayer and in a membrane. 
Note that the initial flat state, consists of two 5's and two 7's, and can be 
regarded as a tightly bound dislocation pair. From Ref. [38]. 



The buckling of a vacancy is illustrated in Fig. 10 p8fl . Even in flat space, the vacancy 
does not have the high symmetry implied by the isotropic impurity source term in Eq. 
( |3.40| ). For the harmonic spring nearest neighbor interaction potential used in Ref. [08| 



an 



initially six-fold symmetric vacancy is "crushed" into the lower symmetry object indicated 
in the figure. The crushed vacancy is equivalent to a tightly bound dislocation pair. Such 
vacancies buckle whenever KqOq/ 'k > 26 f38fl . Figure 10b shows a buckled vacancy for 
K al/ k = 92. Note that the 5-fold disclinations have buckled on opposite sides of the 
membrane. 
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It is also interesting to consider buckling of grain boundaries |[38|| . Three-dimensional 
crystals often consist of randomly oriented grains separated by defect walls. Similar poly- 



crystalline order may appear in partially polymerized membrane vesicles [33,^5,53]. In two 
dimensions, low-angle grain boundaries can be modeled by a row of dislocations, with av- 
erage Burgers vector perpendicular to the boundary Consider first a grain boundary 
in a monolayer. Since all dislocations have the same sign, one might think that the elastic 
energy would be enormous. However, cancellations in the long-range part of the strain field 



insures a finite elastic energy per unit length |55| . Let 9 be the tilt angle relating the mis- 
matched crystallites on either side of the grain. A 9 = 21.8° grain boundary is shown in 
Fig. 11a. The spacing h between dislocations (associated with the large dots in Fig. 11a) is 
given by Frank's law |55|, h = |6sin |. On length scales large compared to h, the composite 
stress field of all the dislocations dies off exponentially, provided the Burgers vectors are 
strictly perpendicular to the boundary; on shorter scales, however, the behavior of the grain 



boundary should be dominated by the indivi 



dual dislocations. 



(a) 



(b) 

FIG. 11. Flat (a) and buckled (b) 21.8° grain boundaries. The two grains 
meet along the row of heavy dots. From Ref. [38]. 

Consider now a grain boundary inserted into a membrane. We have seen that isolated 
dislocations buckle whenever R > Rb ~ 121k/ K^h. One might expect a similar buckling in 
grain boundaries whenever h > Rt>. Numerical studies [33] show that Eq. ( |3.46| ) is again 



satisfied, with £ = ybh and 7 ~ 120, consistent with this guess. A buckled 21.8° grain 
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boundary is shown in Fig. lib for Kobh/n = 300. Note that the membrane remains flat far 
from the boundary. 



D. Defects in Hexatic and Liquid Membranes 

As discussed above, a finite concentration of unbound dislocations is present in mem- 
branes at any nonzero temperature. Provided disclinations remain bound together, the 
resulting phase is a hexatic, similar to the hexatic liquid which arises in the theory of 



two-dimensional melting of monolayers |32|| . Hexatic fluids are characterized by extended 
correlations in the complex order parameter ipeif = exp[6i#(r)], where 9(r) is the angle a 
bond between neighboring particles centered at r makes with respect to a local reference 
axis. Because this reference axis changes when parallel transported on a curved surface, 
there is an important coupling between the bond angle field 9(r) and the geometry of mem- 
brane. This coupling stiffens hexatic membranes relative to their isotropic liquid membrane 
counterparts, although they still fluctuate more wildly than a crystalline tethered surface. 
The continuum elastic-free energy for hexatic membranes was derived in Ref. [35|: 

Fh = \K A J d 2 r - -es&MMdjf))] 2 + i« J d 2 r(V 2 f) 2 (3.47) 

where K A is the hexatic stiffness constant controlling fluctuations in 9(f). Note the similar- 
ity between Fh and the crystalline membrane energy Eq. ( |3.9|) . Both free energies contain a 
bending rigidity, and the out-of-plane displacements act like a "vector potential" coupled to 
the low energy Goldstone modes associated with the relevant broken symmetries. Hexatic 
membranes have many fascinating properties, including a "crinkled" phase, and an interest- 
ing literature is developing around them [56-58] . Here we discuss the von Karman equations 
for hexatics and show that they predict a buckling transition for positive disclination defect 
at a critical value of the ratio n/ K A [39,59, 60]. We argue that the energy for negative discli- 
nations must be different, and point out that this asymmetry has important consequences 
for liquid membranes once disclinations proliferate. 

The variation of Fh with respect to / and 9 leads to 

kV 4 / = K A (d i d J f)t jk d k (d i 9 - A,) , (3.48) 

and 

di(di9 - A) = 0, (3.49) 
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with 



A = -e ]k d k ((dj)(d,f)) . 



(3.50) 



In analogy to Eq. ( |3.16| ), we introduce a hexatic stress function xh (the Cauchy conjugate 
function for the hexatic "current") via 



K A {d,fi - Ai) = eijdjXH 



(3.51) 



so that Eq. ( |3.49| ) is satisfied automatically. Now apply the operator e^dk to Eq. ( |3.51| ) and 
use Eq. (|3.22j ) to rewrite the derivatives of 9 in terms of the disclination density S(r). The 
resulting equation, when combined with Eq. (3.48), results in the "von Karman equations 
for hexatics," 131 



„4, d 2 XH d 2 f d\ H d 2 f d 2 XH d 2 f 

KV J — 1 2 

dy 2 dx 2 dx 2 dy 2 dxdy dxdy 



(3.52a) 



&fd 2 f | / &f\ 2 | 
dx 2 dy 2 \dxdy J 



(3.52b) 



Note that these equations are identical to Eqs. (3.38) except that Kq — > Ka and V 4 - 
in Eq. (3.52b). 

Suppose S{r) = (+ir/3)5(r), representing a single positive disclination at the origin. We 
then look for functions Xui^) and f(f) of the form 



Xh{t) = -«ln(r/a ), f(r) = ±ar 



(3.53) 



where a allows for conical dislocation buckling and remains to be determined. It is easy to 
check that Eq. (|3.52a| ) is obeyed for arbitrary a. All three terms in Eq. ( |3.52b| ) are now 
proportional to delta functions, and equating the coefficients determines a, 

* 2 = 5 " £ ' < 3 ' 54 > 

For k/Ka > 1/6, there are no real solutions and the membrane remains flat. For k/Ka < 
1/6, however, the disclination buckles [39[|. Using Eq. (|3.47|) we readily find that the energy 
of this hexatic disclination in a membrane of radius R is 

|7tk(1 — 3k/Ka) hi(-R/ao)) k/K a < |, 
Z§L\n(R/a Q ), k/K a >\. 



E R 



(3.55) 
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When k/Ka > 1/6, we recover the elastic energy for disclinations in flat hexatic monolayers 
32fl . When Ka — > oo, we recover the buckled disclination energy ( 3.42|) for crystalline 



membranes in the inextensional limit. 
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FIG. 12. Locus of entropic instabilities for a 5-fold disclination in a 
hexatic membrane (solid line). The shaded hexatic region is divided into 
regions where 5-fold disclinations buckle or remain flat, depending on the ratio 
Ka/k. The dashed line is an estimate of where an isolated 7-fold disclination 
becomes unstable. The region above and to the right of the solid line is an 
isotropic fluid. 

The total free energy of the 5-fold disclination including the positional entropy is 



F 5 = E 5 - ST 

= E 5 - 2T\n(R/a ) . 



(3.56) 



This free energy becomes negative above the solid line in Fig. 12. A curve of this kind was 
first presented by Guitter and Kardar [57], who used a more accurate theory to discuss the 
consequences of 5-fold disclination buckling. Their criteria for buckling and for disclination 
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proliferation are similar, but not identical to ours. Small errors are to be expected in our 
approach because the Gaussian curvature G(r) in the Monge representation is actually 

= mmn- wf f (3 . 57) 

(1 + |V/| 2 )V2 

in contrast to the small |V/| approximation which appears in Eq. (3.52b) [60]. A factor (1 + 
|Yj|2^i/2 s } lou i c [ a i so De included in the measure [60]. Nevertheless, the essential predictions 
of the two approaches are identical: the hexatic becomes unstable to an entropically driven 
proliferation of 5-fold disclinations whenever k/T or Ka/T become sufficiently small. 

What about negative disclinations? When Ka — > oo the similar approximations lead 
a result identical to the inextensional crystalline 7-fold defect described by Eq. ( p.44 ). 
Although the energies of both 5- and 7-fold defects diverge logarithmically, the coefficient for 
negative disclinations is larger. This asymmetry persists (but is reduced slightly) in more 
accurate numerical computations [[37]]. When k —>■ oo, we approach the monolayer limit, and 
the 5- and 7-fold disclination energies will diverge logarithmically with equal coefficients. It 
seems plausible that the energy of a 7-fold disclination diverges logarithmically with system 
size for arbitrary values of k/Ka, just as for the 5-fold disclinations. With the two above 
limits in mind, it seems clear that the free energy F 7 = E 7 — 2T\n(R/ao) for negative 
disclinations only becomes negative above the dashed curve in Fig. 12. The energy/entropy 
argument seems to predict that positive and negative disclinations will unbind at two distinct 
temperatures! This conclusion, however, is probably incorrect. Consider the region between 
the solid and dashed lines in Fig. 12. Since the 5-fold defects are energetically favorable, the 
positive-free energy cost of a 7-fold defect can be compensated by a number of nearby 5-fold 
disclinations: A "composite defect" consisting of, say, one 7 and three 5's will have a negative 
total free energy. Thus both positive and negative disclinations will proliferate above the 
solid line. A more complete theory of the disclination transition would include interactions 
between 5's and 7's and might also require overall disclination "charge neutrality." Progress 
in this direction has recently been made by Park and Lubensky [58], although an explicit 
disclination asymmetry has not yet been incorporated into the theory. 

The fundamental asymmetry between positive and negative disclinations should persist 
when hexatics melt into isotropic liquid membranes. Let us denote the defect core energies 
by E§ and E-j and assume that all long-range elastic energies have been screened out. Assume 
further that there is no constraint of disclination charge neutrality. This would be the case in 
membranes with free edges, or for the experimentally relevant case of liquid bilayer surfaces 
which can change their genus freely. The areal densities of disclinations will then be different, 
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n 5 « Oq 2 e" jB5/T , n 7 « a^ 2 e~ E7/T . (3.58) 

Which defect predominates in liquid membranes depends on microscopic details such as 
interaction potentials, etc. Close to a transition to a hexatic phase, however, the arguments 
given above suggest that 5-fold disclinations dominate. 

Disclination asymmetry in liquid membranes parallels the behavior of vacancies and 
interstitials in conventional crystalline solids p8| . Although vacancies and interstitials are 



the antidefects of each other, they nevertheless have very different energies. Unless vacancies 
and interstitials are created from a perfect crystal with periodic boundary conditions, they 
will in general occur with different concentrations. Periodic boundary conditions create an 
artificial constraint which forces vacancies and interstitials to be created in pairs rather 
than diffusing in from the surface. The concentrations must also be equal if the defects 
are charged, as in ionic crystals ||28|| . The flatness of monolayers similarly constrains the 



disclination densities to be equal at all stages in the theory of two-dimensional melting |32 
In membranes artificially constrained to have the topology of a torroidal surface, the numbers 
of 5- and 7-fold disclinations are forced to be equal to Euler's theorem. A calculation similar 
that for point defects in ionic crystals then leads to 

n 5 = n 7 « afe-^ +E ^l 2T . (3.59) 

The free energy of liquid membranes is often expressed in the Helfrich form [62], 
1 



. , 2 k J d 2 r(V 2 fY + k g J dV[(fl£/)(fig/) - (<W)1 • (3-60) 

The first term is the usual bending rigidity. The remaining one is proportional to the 
integrated Gaussian curvature and its coefficient kg is often called the Gaussian rigidity. 
This second term is a perfect derivative, and integrates to a constant for surfaces of fixed 



genus pi |. The microscopic origins of kq are obscure. Its sign, however, clearly determined 
by the disclination asymmetry discussed above: An excess of 5-fold disclinations corresponds 
to kg < and favors positive net Gaussian curvature. Membrane phases many spherical 
vesicles will predominate in this case. An excess of 7-fold disclinations means > 0, and a 
bias toward negative Gaussian curvatures. Complex "plumbers nightmare" lipid membrane 
phases [34] are then favored. 

On a more formal level, it must be the case from Euler's theorem that the Gaussian 
curvature integrated over a membrane with free edges of area Q gives the disclination asym- 
metry 

Jd 2 r^gG(f)=Q(n 5 -n 7 ), (3.61) 
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where Jg = 1 + |V/| 2 and G{r) is given by Eq. (|3.57| ). Equation (|3.60|) can thus be rewritten 



as 



where N 5 and N 7 are the total numbers of 5- and 7-fold disclinations in the membrane. The 
Gaussian rigidity Kg thus acts as a chemical potential which must be adjusted to give the 
correct asymmetry between the populations of 5's and 7's. In this sense, its effect is similar 
to that of a nonzero magnetic field in Eq. (|2.56|) , which would lead to a net excess of positive 
or negative vortices in superconducting films. For a related perspective on the physics of 
metallic glasses, see Ref. p2f . 

It would be interesting to search for the disclination asymmetry discussed here in com- 
puter simulations of membranes. The most straightforward approach would be to study 
initially flat (i.e., large k) liquid membranes with free edges, so that disclinations can enter 
and exit freely at the boundary. When the bending rigidity is reduced, a bias in the aver- 
age Gaussian curvature should emerge as the membranes curl up into the third dimension. 
One could then vary the interparticle potentials to study what factors influence the sign of 
Ei — E§. Such an understanding might lead to controlled synthesis of membranes with a 
predetermined sign of k. 
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APPENDIX: SUPERFLUID DENSITY AND MOMENTUM CORRELATIONS 
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FIG. 13. Two experiments which lead to the superfluid density. In (a), 
there are periodic boundary conditions in the x direction, and impenetrable 
walls at y = ±L y /2. These boundary conditions are reversed in (b). In both 
cases, the walls and substrate move at velocity u x in the x direction. 

Consider a simplified version of the oscillating substrate experiment discussed in Sec. IIA. 
Imagine that the substrate is wrapped around to form a cylinder, coated uniformly with 
He 4 . We neglect substrate inhomogeneities, as is appropriate near the transition or if the 
films are sufficiently thick [63]. Now imagine that the substrate oscillations are around the 
cylinder axis, and slow enough so that the motion of the substrate is essentially a uniform 
translation with velocity u. As illustrated in Fig. 13a, we choose a coordinate system such 
that the substrate moves along the x axis. The cylinder has height L y along the y axis, 
and there are barriers at y — ±L y /2 which prevent the film from escaping in this direction. 
Its circumference is L x . In two or three dimensions, the superfluid density measures the 



response of liquids to moving walls [20-22]- In Fig. 13a, the "walls" are provided by the 
substrate itself plus the barriers at the top and bottom of the cylinder. Instead of solving 
a complicated statistical mechanics problem with moving boundaries, it is easier to make a 
Galilean transformation to a coordinate system which moves at the velocity of the substrate, 
so that the boundaries are fixed. In the original laboratory frame of reference, the superfluid 
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fraction of the film will remain at rest. Hence, it appears to be moving with average velocity 
— u after the Galilean transformation. 

If Ji is the Hamiltonian in laboratory frame, averages after the Galilean transformation 
must be computed with respect to the Hamiltonian [64]. 

-H' = n-H. p + - Mu 2 , (Al) 

where P is the total momentum of the film and M is the total mass. Upon coarse graining 
the film, the free energy which appears in Eq. ( |2.10| ) is correspondingly replaced by 

F' = F-u - [ g{r)d 2 r + 0{u 2 ) , (A2) 

where g(r) is given by Eq. (2.43). The relative weights of different configurations of ip(r) in 
the moving frame with stationary boundary conditions, given by exp[— F'/T], must be the 
same as those in the laboratory frame with moving boundary conditions. The average value 
of the momentum in the laboratory frame is thus given by 

where g = (h/2i)[ip*(r)V'ip(r) — ijj(r) Vip{r)\ and the subscript on the average denotes a 
system with boundaries moving at velocity u. To linear order in the wall velocity, we then 
have 

( 9i (f)U = f I d 2 r' (gm 9j (f '))n=o Uj + 0{u 2 ) , (A4) 

and we see that the momemtum generated by moving walls in a helium film is determined 
by the momentum correlations for a system with the walls at rest. We now expand g(r) in 
Fourier variables, 

m = ^T,9(q)e^, (A5) 

with 

g(q) = J d 2 r e^g^ , (A6) 
and where fi is the system area. Equation (A4) can now be rewritten as 



(<fc(r))a = \ hm a^qjUj + 0{u 2 ) (A7) 

1 q^O 



where 
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Cij(q) = (gi(q)9j(-q))u=o 

= A( q )^ + B(q)(6 v -^ , (A8) 

and A(q) and B(q) are functions only of the magnitude of q. 

The coefficient of Uj Eq. ( |A7| ) involves a delicate limiting procedure whenever A{q) ^ 



B(q) p0|-|2^|. Consider first the situation outlined above, with the periodic boundary condi- 
tions in the x-direction appropriate to a smooth cylindrical substrate. For any finite cylinder 
circumference L x , there are always Fourier components of the g(r) at q x = 0, projected out 
by the x' integration in Eq. (|A4j) . Along the y direction, however, the q y = mode only 
appears in the limit L y — > oo, because g(r) must vanish at the top and bottom of the cylin- 
der. Thus, the correct order of limits for large sample sizes in this experiment is the limit 
q x — > 0, followed by the limit q y — > 0. The x component of the film momentum induced by 
the moving substrate is evidently 

(9x{r))u = ~ hm hm C xx (q)u x 

= ^«.. (A9) 

Since only the normal component of the film moves with the substrate, we identify p n (T) 
with the coefficient of u x , 

Pn(T) = ^ . (A10) 

Now consider a related, but fundamentally different experiment: We repeat the cylinder 
periodically along the y axis, but erect impenetrable barriers along the x axis at fixed 
positions x = ±L x /2. See Fig. 13b. The new barrier (equivalent to gouging a slit in the 
substrate parallel to the axis of the cylinder) prevents the film from circulating completely 
around the cylinder. The correct limiting procedure is now the limit q y — * 0, reflecting 
the periodic boundary conditions along y, followed by the limit q x — > 0, describing the new 
barriers in the limit L x — > oo, 

(9x{r)) = ^ hm hm C xx (q)u x 

1 qx^Q q v — >0 

= fU. (All) 

In this experiment, all of the film must clearly move at the substrate velocity since it is pushed 
around by the barrier. Now we identify p to t, the total film density, with the coefficient of 
u x , so that 
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P» = ^ . (A12) 

In a normal liquid, these two different limiting procedures would lead to identical results. 
In a superfluid, however, there is a nonzero superfluid density, defined by 

Ps(T)=p tot -p n (T) . (A13) 

The formula for the superfluid density is thus 

p s {T) = 1 lim [A(q) - B(q)\ . (A14) 
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